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SUMMARY 


This  report  contains  an  analysis  of  the  experimental  results  obtained 
from  four  supercritical  natural  laminar  flow  (NLF)  airfoils  investigated 
in  the  NAE  High  Reynolds  Number  Test  Facility.  The  airfoils  have 
maximum  thickness  to  chord  ratios  of  0.1 0.13,  0.16  and  0.21  and  were 
designed  for  a  lift  coefficient  of  Cl  =  0.6.  .eir  design  Mach  numbers  were 
0.8,  0.76,  0.72  and  0.68  respectively  and  the  design  chord  Reynolds  number 
was  12.5  million.  It  was  found  that  all  the  airfoils  showed  the  presence  of  a 
drag  bucket  close  to  design  conditions  and  long  ler  -ths  (in  some  cases 
about  70%)  of  natural  laminar  flow  at  Reynolds  oer  6.7  million.  The 
minimum  drag  for  the  airfoils  was  found  to  u  nge  i  jm  0.0045  to  0.0051, 
representing  far  lower  levels  than  any  airfoil  dominated  by  turbulent 
boundary  layer.  It  is  also  indicated  that  with  transition  fixed  at  about  10% 
chord  the  drag  levels  were  similar  to  other  airfoils  with  turbulent  boundary 
layers. 


RESUME 

Le  present  rapport  contient  une  analyse  des  resultats  experimentaux 
obtenus  pour  quatre  profils  a6rodynamiques  h  ecoulements  laminaires 
naturels  surcritiques  qui  ont  fait  l'objet  de  recherches  h  l'installation 
d'essai  pour  les  nombres  de  Reynolds  61ev6s.  Les  profils  ont  des  rapports 
epaisseur/corde  maximaux  de  0,10,  0,13,  0,16  et  0,21,  et  sent  con?us  pour  un 
coefficient  de  portance  de  Cl  =  0,6.  Leurs  nombres  de  Mach  theoriques  sont 
de  0,8,  0,76,  0,72  et  0,68  respectivement,  et  le  nombre  de  Reynolds  theorique 
des  cordes  est  de  12,5  millions.  Les  resultats  rev&lent  que  tous  les  profils 
presentent  un  godet  de  trainee  au  voisinage  des  conditions  theoriques  et  des 
longueurs  importantes  (parfois  d'environ  70  %)  d'6coulement  laminaire 
naturel  au  nombre  de  Reynolds  de  6,7  millions.  La  trainee  minimale  pour 
les  profils  varie  entre  0,0045  et  0,0051,  ce  qui  est  de  beaucoup  inferieur  5  celle 
de  tout  profil  domine  par  une  couche  limite  turbulente.  A  noter  aussi  que, 
pour  une  transition  fixee  a  environ  10  %  de  la  corde,  la  trainee  est  voisine 
de  celle  d’autres  profils  avec  couches  limites  turbulentes. 
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1 .  INTRODUCTION 

Future  transport  aircraft  will  benefit  from  improved 
airfoil  designs  that  reduce  wing  section  drag.  Toward  that 
objective,  the  National  Aeronautical  Establishment  and  Boeing 
Canada  de  Havilland  Division  have  had  an  ongoing  program  of 
research  and  development  aimed  at  developing  improved 
supercritical  airfoils  suitable  for  future  regional  transport 
aircraft.  The  purpose  of  this  report  is  to  provide  a  compendium 
of  the  results  obtained  over  several  years  of  designing  and 
testing  four  supercritical  airfoils  that  are  capable  of 
supporting  extensive  regions  of  natural  laminar  flow  { NLF )  in 
suitably  favorable  conditions.  This  technology  will  provide 
substantial  reduction  in  drag  or  alternatively  allow 
substantial  increases  in  wing  thickness  for  a  given  Mach  number 
and  lift  coefficient. 


when  designing  these  airfoils,  one  of  the  foremost 
objectives  was  to  ensure  that  good  aerodynamic  characteristics 
were  retained  when  the  boundary  layers  were  made  turbulent  from 
near  the  noses  of  the  airfoils.  This  situation  can  be  expected 
to  prevail  occasionally,  for  premature  transition  can  be  caused 
by  some  type  of  surface  contamination  or  free  stream 
turbulence.  However,  for  circumstances  where  conditions  are 
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are  designed  to  support  laminar  flow  back  to  45  to  70%  chord  on 
their  upper  surfaces  and  up  to  40  to  50%  chord  on  their  lower 
surfaces.  This  extent  of  laminar  flow  will  lead  to  drag 


reduction  of  at  least  50%  relative  to  the  same  airfoils  with 
turbulent  boundary  layers. 


Validating  the  performance  of  NLF-capable  airfoils  at 
high  Reynolds  number  becomes  difficult  in  wind  tunnels  due  to 
the  need  for  very  low  turbulence  levels  and  for  extremely  fine 
tolerances  on  small  models.  Measurements  made  in  1978  by 
Elf strom  in  the  NAE  Two  Dimensional  Test  Facility  [1]  showed 
turbulence  levels  were  about  0.30%  at  Reynolds  numbers  of  14-20 
million  .  Calculations  made  using  this  level  of  turbulence 
with  the  method  of  Van  Driest  and  Blumer  [2]  suggested  there 
was  little  prospect  for  extensive  NLF  under  these  conditions. 
However,  subsequent  tests  at  NAE  showed  lower  turbulence  levels 
of  about  0.16%  as  a  result  of  new  screens  being  installed  [3] 
and  edgetone  noise  effects  from  the  porous  walls  being  reduced. 


Four  airfoils  were  designed  and  tested.  They  ranged 
from  10%  to  21%  in  maximum  thickness  to  chord.  Their  geometries 
are  shown  in  Fig.l.  The  airfoil  designations  sho.vn  in  Fig.l 
e.g.,  NAE  76-060-13(1)  indicate  the  design  Mach  number,  design 
lift,  maximum  thickness  to  chord  and,  in  brackets  the  version 
of  the  airfoil.  The  design  lift  coefficient  selected  for  the 
airfoils  was  0.6  to  suit  aircraft  operating  over  long  ranges 
and  at  high  altitudes.  The  airfoils  were  designed  using  BGK 
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and  accelerated  flow  up  to  the  shock  wave,  in  anticipation  that 
accelerated  flow  would  lead  to  delayed  transition  and  extended 
runs  of  laminar  flow.  Note  that  the  trailing  edge  thickness 
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(~0.5%)  of  the  16  and  21%  airfoils  is  larger  than  that  of  the 
thinner  air  foils ( ~0 . 1%  )  .  The  thinner  trailing  edge  for  later 
designs  was  believed  to  be  more  advantageous  for  NLF  design  as 
observed  in  Reference  4. 

Some  features  of  the  design  velocity  distribution  for 
the  airfoils  are  summarized  in  Fig.  2.  At  the  design  condition 
the  peak  Mach  numbers  on  the  upper  surface  were  kept  to 
1.10-1.12  to  give  low  wave  drag,  and  the  velocity  gradients 
were  made  slightly  favorable  to  encourage  NLF.  A  moderate 
amount  of  aft  loading  was  used  to  enhance  performance,  subject 
to  the  constraints  that  pitching  moment  coefficients  should  not 
be  less  than  -0.14  at  the  design  conditions.  Boundary-layer 
calculations  were  made  with  almost  fully  turbulent  conditions 
which  showed  good  margins  from  flow  separation  on  the  upper  and 
lower  surfaces  of  the  airfoils  at  their  respective  design 
conditions . 

The  airfoils  were  tested  in  the  NAE  0.38m  X  1.5m 
two-dimensional  wind  tunnel  over  a  range  of  chord  Reynolds 
numbers  from  6.7  to  16.7  million  .  The  tunnel  is  a  blowdown 
type  with  the  top  and  bottom  walls  having  normal  holes  giving  a 
porosity  of  20.5%,  and  the  sidewalls  having  controlled  suction 
in  the  region  of  the  model.  The  model  of  the  16%  t/c  (thickness 
to  chord )  airfoil  v/  a  s  made  of  12  inches  chord  and  f  h  a  models  of 
the  other  airfoils  were  10  inches  chord.  Fig.  3  shows  the 
schematics  of  the  2  D  test  section  (0.38  m  X  1.5  m)  which  is 
mounted  in  the  1.5  m  X  1.5  m  test  section  of  the  NAE  trisonic 
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wind  tunnel  for  these  tests.  See  Reference  5  for  more  detailed 
description  of  this  facility. 

It  should  be  noted  that  all  airfoils  were  tested  at 
the  same  unit  Reynolds  number.  Thus  the  chord  Reynolds  number 
for  the  16%  thick  airfoil  is  1.2  times  larger  than  that  of  the 
other  three  airfoils.  In  the  following  text,  the  Reynolds 
numbers  quoted  are  referred  to  the  10  inch  chord  airfoils, 
except  otherwise  indicated,  the  Reynolds  number  for  the  16% 
thick  airfoil  is  1.2  times  the  cited  value. 

This  report  will  first  review  the  experimentally 
determined  characteristics  of  all  the  airfoils  under  both  fixed 
transition  (usually  7%  and  15%  on  the  upper  and  lower  surfaces 
respectively)  and  free  transition  conditions.  A  later  section 
will  compare  computed  results  from  a  recent  code  GRUMFOIL  with 
the  experimental  data  for  the  10  and  13%  airfoils. 
Unfortunately,  transition  free  computed  results  on  the  13% 
airfoil  are  not  available  as  convergence  problems  were 
encountered  with  the  GRUMFOIL  computer  code.  Other  comparisons 
of  theory  with  experiment  for  all  the  airfoils  under  transition 
fixed  conditions  can  be  found  in  Ref.  6  while  Ref.  7  gives  more 
information  on  the  two  thicker  airfoils. 

2.  DT  SCtlR  ST  ON  OF  KXPERTMF.NTAT.  RESIIT/PS 
2.1  LIFT  PERFORMANCE 


2.1.1  Lift  versus  angle  of  attack 
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Figure  4  shows  the  transition  free  CLB  versus  a 

performances  for  all  the  four  airfoils  at  Mach  number  of  M=  0.3 

and  Reynolds  number  of  6.7  X  10^  (  16%  thick  airfoil  at  M=  0.2 

6  —  1 

and  Re/c  =5.9  X  10  ft  ).  CL  MAX  appears  to  increase  with 
the  thickness  of  the  airfoil  starting  from  a  low  of  about  1.19 
at  a=  9.18  deg  for  the  10  %  thick  (t/c)  airfoil,  to  the  highest 
CT  value  of  1.68  at  a  =  13.67  deg  for  the  21  %  thickness 
airfoil.  The  Cr  value  for  the  10%  airfoil  is  well  below  the 

other  thicker  airfoils.  The  leading  edge  of  the  10%  airfoil  has 
a  nose  radius  much  smaller  than  the  other  thicker  airfoils,  and 


as  explained  in  Ref.  [8],  a  decrease  in  the  nose  radius  can  be 


responsible  for  a  substantial  drop  in  the  C 


L  MAX 


value  for  the 


same  airfoil. 


In  Fig.  4  the  CL  MAX  values  in  most  cases  are  well 
defined  and  the  C. n  versus  a  are  fairly  linear.  At  a  =  0.0  the 
CLB  values  range  from  0.275  for  10%  thick  airfoil  to  0.445  for 
21%  thick  airfoil.  The  increase  must  result  from  increasing 
camber  for  thicker  airfoils,  see  Fig  1.  At  other  Reynolds 
numbers,  the  angle  of  attack  was  not  traversed  up  to 
sufficiently  high  a's  to  provide  a  complete  C.  R  v  a  picture. 


2.1.2  Lift  curve  slope  versus  Mach  number 


Figures  5  and  6  show  the  transition  free 

cases  of  dC^/da  (at  CL  =  0.0)  against  Mach  number  for  the  four 

6  6 

airfoils  at  Reynolds  numbers  of  6.7  X10  and  12.5  X  10 
respectively.  From  Figure  5  it  appears  that  dC^/da  values  remain 
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somewhat  steady  up  to  M  =  0.5  after  which  point  they  show  an 
increasing  trend,  reaching  a  maximum  at  or  close  to  the 
individual  design  Mach  number  of  the  airfoil.  This  increasing 
trend  with  Mach  number  is  consistent  with  Prandtl-Glauer t 
theory  (compressibility  effect)  which  predicts  a  linear 
increase  with  respect  to  the  parameter  =  1//(1-M  )  when  close  to 
transonic  Mach  Numbers.  Similar  trends  have  been  observed  for  a 
NACA  230  series  airfoil  tested  in  DVL  (Deutsche  Versuchsanstalt 
fur  Luftfahrt)  2.7m  wind  tunnel  Ref  [8].  The  3CT/3a  data  for  the 

g 

higher  Reynolds  number  case  12.5  X  10  are  somewhat  lower  than 
those  at  the  lower  Reynolds  number  6.7  X  10^  .  The  thicker 
turbulent  boundary  layer  effectively  decambers  the  airfoil 
leading  to  lower  3CL/Sa. 

To  investigate  linear  compressibility  similitude 

amongst  the  four  airfoil  ,  3CL/3a  data  at  Re  =  6.7  and  12.5 

million  under  free  transition  conditions  were  also  plotted 

against  1//(1-M  ).The  data  shown  in  Fig  7  seem  to  collapse  well 

for  M  <  0.64  suggesting  similarity  in  performance  at  subsonic 

Mach  numbers.  Note  also  that  3CL/3a  is  fairly  linear  with  the 

2 

parameter  1//(1-M  )  up  to  an  abscissa  value  of  about  1.25. 

Figures  8  and  9  are  the  plots  for  the  transition 
fixed  values  of  3CL/3a  against  Mach  number  at  Reynolds  numbers 
of  6.7  X  10^  and  12.5  X  10^  respectively.  There  is  no  data 
available  for  M  <  0.6.  For  the  Mach  number  range  covered  in 
these  plots  0.6  <  M  <  0.8  ,  in  general  there  is  an  increase  in 
3CL/3a  values  with  Mach  number,  with  some  airfoils  showing  a 
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distinct  maximum  towards  the  upper  end  of  the  range.  Under 

transition  fixed  conditions  it  is  noted  that  for  a  given  Mach 

number  the  thinner  boundary  layer  at  Re/c  =  12.5  XlO^  gives 

rise  to  higher  3CL/3cc  values  than  their  counterparts  at  6.7  X 

10b  .  Linear  compressibility  was  again  studied  by  plotting 

3Cr/3a  against  1//{1-M  )  for  the  transition  fixed  data  in  Pig 
w 

10.  Again,  the  data  seem  to  collapse  only  for  the  subsonic  Mach 
numbers  M  <  0.67. 

It  should  be  mentioned  that  von-Xarman  [3]  transonic 

similitude  was  also  investigated  for  the  above  data  under 

transition  fixed  and  free  conditions.  This  exercise  was  carried 

?  ')  n  ^  3 

cut  by  plotting  3CL/3<x(  (  y+1  )M  /S  )  against 

( 1-ri2 )/(  S(  y+1  )M2 ) 0  ‘ 66^  (  For  the  sake  of  brevity  this  plot  is 

not  shown).  No  satisfactory  collapse  of  data  for  the  four 
airfoils  was  noted. 

Comparison  of  transition  free  data  at  Re/c  =  6.7  X10® 
(Figure  5)  against  the  corresponding  turbulent  airfoil  data 
(Figure  8)  shows  that  there  is  a  decrease  in  3CL/3a  values  once 
the  laminar  flow  is  lost.  However,  at  the  higher  Reynolds  number 
of  12.5  X  10  ,  the  difference  in  (compare  Figures  6  and  9) 

3CT/3a  values  for  the  transition  free  and  transition  fixed  runs 

li 

is  not  as  noticeable. 

2.1.3  Maximum  lift  and  lift  at  separation  onset 
versus  Mach  number 
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Figures  11  and  12  show  the  free  transition 
lift  performance  of  the  four  airfoils  against  Mach  number  at 
chord  Reynolds  numbers  of  6.7  X  10^  and  12.5  X  10^ 
respectively.  For  each  airfoil  the  C.  performance  is  also 

L  SPp 

shown  along  with  the  CT  .  C.  in  each  case  is  obtained  by 

Li  max  Li  sap 

plotting  C  at  or  close  to  a  chord  station  x/c  =  0.96  versus 

r 

CrD,  and  then  determining  the  point  on  the  curve  where  3C  /8CT  « 

Lb  p  L 

-0.4,  Ref. 10.  Lift  corresponding  to  the  drag  rise  Mach  numbers 
was  also  evaluated  (  based  on  a  3CD/9M  =0.1  method)  from 
appropriate  drag  polars  and  is  also  shown  in  Figures  11  and  12. 

When  discussing  the  performance  of  these  airfoils  as  a 

group  it  must  be  understood  that  each  one  of  them  has  been 

designed  to  perform  best  at  its  own  design  conditions.  For 

example  at  M  =0.7  the  21%  airfoil  is  at  or  close  to  its  design 

lift  condition  whereas  the  10  %  airfoil  is  far  from  its  design 

Mach  number  value.  Equally  the  10%  is  designed  for  M=0.8,  where 

its  lift  performance  ij^,  and  expected  to  be  better  than  the 

other  three.  In  general  a  glance  at  Figure  11  indicates  that  for 

all  the  airfoils  both  Cr _ and  CT  „  decrease  with  increasing 

L  max  L  sep  3 

Mach  number.  Another  interesting  point  to  note  is  that  at  higher 
Mach  numbers  M  >  0.7  ,  there  is  a  tendency  for  thinner  airfoils 
to  show  higher  CT  „„  values  and  vice  versa  for  M  <  0.7  which  is 
indicative  of  different  lift  developing  mechanism  in  different 
Mach  number  regimes. 

The  drag  rise  Mach  number  data  in  both  Reynolds  number 
cases  (Figures  11  and  12)  also  conforms  to  the  design 
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constraints  of  each  airfoil  and  cannot  be  discussed  in  a 
comparative  sense. 

Figures  13  and  14  show  the  corresponding  transition 
fixed  data  for  the  four  airfoils  at  chord  Reynolds  numbers  of 

ft 

6.1  X  10  and  12.5  X  10  respectively.  The  general  behavior  of 

CT  „  and  CT  again  shows  a  decreasing  trend  with  the 

L  max  L  sep  3  3 

increase  of  Mach  number  with  the  overall  levels  of  lift  (in  most 
cases)  being  a  little  bit  below  the  transition  free 
counterparts.  There  is  no  radical  shift  in  drag  rise  Mach  number 
in  going  from  transition  free  to  turbulent  conditions  for  the 
two  thicker  airfoils,  however,  the  thinner  13  %  and  10  % 
airfoils  do  appear  to  show  a  drop  in  drag  rise  Mach  number  when 
the  turbulent  flow  conditions  are  forced. 


As  expected  in  almost  all  of  the  above 

Mach  number  plots  CL  sgp  values  always  trail  the 

cr  data. 

L  max 


lift  versus 
corresponding 


2.2  PITCHING  MOMENT  RESPONSE 


2.2.1  Pitching  Moment  versus  Mach  Number 

Figures  15  and  16,  respectively  show  the 
transition  free  and  transition  fixed  pitching  moment  versus  Mach 

number  of  the  four  airfoils  at  CLg  =  0.6  and  chord  Reynolds 

£ 

Number  of  6.7  X  10  .  For  the  transition  free  results  in  Figure 
15  ,  it  is  clear  that  both  the  thickness  (t/c)  and  the  Mach 
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number  have  the  effect  of  increasing  the  magnitude  of  the 
pitching  moment.  As  mentioned  earlier,  at  the  design  Mach  number 
these  airfoils  were  supposed  to  have  a  pitching  moment  greater 
than  -0.14.  With  the  exception  of  the  10%  airfoil  which  appears 
to  have  a  pitching  moment  value  of  about  -0.14  at  M=0.8,  all 
others  fall  well  below  this  constraint  at  their  respective 
design  r  ch  numbers.  The  16  %  airfoil  (NAE  72-060-16:1)  perhaps 
displays  the  largest  pitching  moment  magnitude  | |  in  almost 
the  entire  Mach  number  range  (  with  the  exception  of  a  small 
interval  0.68  <  M  <  0.72  in  which  the  21%  airfoil  takes  over). 


The  transition  moves  forward  with  the  increase  of  the 
Reynolds  number  and  a  thicker  boundary  layer  with  a 
comparatively  smaller  loading  at  the  aft  portion  causes  a 
reduction  in  the  pitching  moment  magnitude  as  shown  by  the 
comparison  of  the  turbulent  airfoil  data  in  Figure  16  to  the 
free  transition  data  in  Figure  15.  The  pitching  moment 
magnitude  still  shows  a  similar  trend  with  thickness  and  Mach 
number  as  was  observed  for  the  transition  free  case  in  Figure 
15. 


Figure  17  and  18  respectively  show  a  similar  transition 
free  and  transition  fixed  comparison  at  CTD  =  0.6  and  a  still 

LB 

c. 

higher  chord  Reynolds  number  of  12.5  X10  .  The  transition  free 
data  in  Figure  17  once  again  show  that  the  airfoils  thickness 
and  increase  of  Mach  number  have  the  effect  of  increasing  the 
magnitude  of  the  pitching  moment.  Note  that  the  16  %  airfoil 
even  at  this  higher  Reynolds  number  of  12.5  X  10  seems  to 


provide  the  largest  pitching  moment  at  almost  all  Mach  number 
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conditions.  It  is  also  interesting  to  note  that  for  the 
transition  free  data  in  Figure  17,  the  thicker  airfoils  21  and 
16  %  appear  to  show  a  bottoming  out  of  the  pitching  moment 
response  followed  by  (in  case  of  21%  airfoil)  a  short  reversal 
caused  by  a  trailing  edge  separation.  Under  transition  fixed 
conditions  in  Figure  18  the  overall  trends  remain  generally 
similar  to  their  transition  free  counterparts. 

2.3  AIRFOIL  DRAG 


The  wake  drag  was  measured  by  the  standard 
sidewall-mounted  traversing  rake  supporting  four  pitot  probes, 
The  method  is  based  on  measuring  the  momentum  deficit  in  the 
wake.  Reference  11  gives  details  of  the  instrumentation.  The 
signal  from  probe  4  is  ignored  as  it  is  sometimes  affected  by 
the  disturbed  sidewall  boundary  layer.  The  total  effective  drag 
CDW  is  usually  computed  from  an  average  of  the  remaining  three 
probes,  see  Figure  (19). 

2.3.1  Drag  Polars 

2. 3. 1.1  Transition  Free  Drag  Polars  at  Re  =  6.7  x  10^ 


Figures  20  to  25  present  transition  free  drag  polars  for 
the  four  si rf oils  st  Re  —  6.7  X  10°  for  the  Ms ch  number  rsn3S 
from  0.3  to  0.7.  For  the  16  %  airfoil  ,  the  low  Mach  number  (  M 
=0.2  )  measurements  were  made  at  Re  =  6  X10^  and  are  thus  duly 
presented  for  comparison. 
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At  M  =  0.3  in  Figure  20  ,  in  the  lower  drag  region  of 
the  polars  with  0.0  <  CLB  <  0.75  ,  the  drag  increases  with  the 
thickness  of  the  airfoil,  with  the  thinnest  10  %  airfoil, 
occupying  the  lowest  drag  level  of  about  0.0046  which  increases 
to  a  value  of  about  0.0080  for  the  21  %  airfoil.  At  M  -0.5  in 
Figure  21,  a  similar  trend  is  repeated  except  in  the  region  0.5 
<  CLB  <  0.75,  where  the  10  %  airfoil  has  very  similar  drag 
characteristics  as  the  13%  airfoil. 

The  drag  polar  for  the  16%  airfoil  was  not  available  at 
M  =  0.6  in  Figure  22.  Drag  polars  for  the  remaining  airfoils  are 
shown.  Except  for  one  point  at  CLQ  =  0.198,  the  drag  polar  for 
the  21%  airfoil  is  notably  higher  than  the  polars  for  the  other 
two  airfoils  which  at  least  in  the  range  0.3  <  CLB  <  0.7  show 
very  similar  drag  levels.  Figure  23  shows  the  drag  polars  for 
the  four  airfoils  at  M  =  0.7.  This  Mach  number  is  very  close  to 
the  design  Mach  number  value  of  the  21  %  airfoil  and  ,  near  the 
design  lift  condition  in  the  range  0.3  <  CTn  <  0.6,  this 
airfoil  displays  the  lowest  drag.  This,  as  will  be  addressed 
later,  is  due  to  extensive  NLF  at  this  condition. 

At  Mach  number  M=0.76  ,  the  drag  polar  for  the  thickest 
21  %  airfoil  was  not  measured  as  the  conditions  here,  are  well 
past  the  drag  rise  Mach  number  for  this  airfoil.  The  polars  for 
the  remaining  three  airfoils,  shown  in  Figure  24,  are  fairly 
similar  and  crisscross  each  other  in  the  range  0.3  <  Crn  <  0.7. 
For  the  thinnest  airfoil  (to  show  the  Mach  number  effect),  the 
drag  polar  at  M  =  0.74  is  also  included  (plotted  using  symbol 
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+).  Only  the  thinnest  two  airfoils  were  tested  at  high  Mach 
numbers  of  0.8  and  0.81.  The  10  %  airfoil  which  at  these  Mach 
numbers  is  performing  at  or  close  to  its  design  conditions  gives 
significantly  less  drag  than  the  13  %,  see  Figure  25.  The  13  % 
airfoil  is  almost  into  drag  rise  but  still  displays  good  drag 
features  for  between  0.2  and  0.5. 

2. 3. 1.2  Transition  Fixed  Drag  Polars  at  Re  =  6.7  X  10^ 

Transition  fixed  drag  polars  are  given  in  Figures  26  to 
29.  In  comparison,  the  data  collected  for  the  transition  fixed 
cases  is  less  detailed  .  The  low  Mach  number  and  Reynolds  number 
case  shown  in  Figure  26  suggests  that  from  the  drag  viewpoint 
there  is  very  little  to  choose  between  the  two  airfoils.  At 
these  flow  conditions  the  boundary  layer  displays  similar 
characteristics  irrespective  of  the  thickness  of  the  airfoil.  At 
a  higher  Mach  number  of  M  =  0.6  in  Figure  27  on  the  other  hand, 
an  increase  in  drag  is  noted  with  the  thicker  airfoils. 


Figures  28  and  29  the  transition  fixed  drag  polars 

at  Mach  numbers  0.68  <  M  <  0.8.  Since  this  Mach  number  range 
includes  the  design  Mach  number  of  all  the  airfoils  ,  the 
minimum  drag  levels  here  are  considerably  higher  than  their 
counterparts  in  transition  free  polars.  In  other  words,  the 
generous  extent  of  natural  laminar  flow  available  is  lost  once 
the  transition  is  fixed  giving  rise  to  the  resulting  wake  drag 
levels.  For  example  the  13  %  airfoil  here  in  Figure  29  at  M  = 


0.76  gives  a  minimum  drag  level  of  0.0102  compared  to  a  minimum 
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drag  level  of  only  about  0.0045  under  natural  laminar  flow 
conditions  in  Figure  24.  It  should  be  stressed  however  that  the 
'turbulent  flow'  drag  levels  are  still  competitive  with  airfoils 
specifically  designed  for  early  transition. 

2. 3. 1.3  Transition  Free  Drag  Polars  at  Re  =  12.5  X  10® 

Figures  30  to  34  show  transition  free  polars  at  a  higher 
Reynolds  number  of  Re  =  12.5  X  10  in  the  Mach  number  range  0.5 
to  0.8.  Since  at  this  high  Reynolds  number  the  airfoils  do  not 
support  as  much  laminar  flow  as  they  did  at  the  lower  chord 
Reynolds  number  case  of  6.7  X  10®  ,  the  airfoils  have 
comparatively  higher  drag  particularly  in  cases  close  to  their 
design  Mach  numbers.  Consider  for  example  the  10  %  airfoil  in 
Figure  34,  which  shows  drag  values  of  about  0.0064  to  0.0080  for 
CLB  <  0.6.  This  drag  level  was  as  low  as  0.0048  -  0.0060  for  the 
same  C^B  limits  at  a  lower  chord  Reynolds  number  of  6.7  X  10® 
in  Figure  25. 

2. 3. 1.4  Transition  Fixed  Drag  Polars  at  Re  =  12.5  X  10® 


Transition  fixed  studies  were  also  done  at  chord  Reynolds 
number  of  Re  =  12.5  X  10  .  The  drag  trends  as  observed  in  the 
drag  polars  for  the  Mach  numbers  M  =  0.6  to  0.8,  shown  in 
Figures  35  to  38,  in  general  remain  the  same  as  under  transition 


free  high  Reynolds  number  polars  in  Figure  30  to  34.  Only  the 
actual  drag  levels  are  somewhat  higher.  Compare  for  instance  the 
transition  fixed  results  for  the  three  thicker  airfoils  at  M  = 
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0.6  in  Figure  35  to  their  transition  free  counterparts  in  Figure 
31.  It  appears  that  even  at  Re/c  =  12.5  X  10^  the  airfoils 
carry  some  laminar  flow  which  is  lost  when  the  transition  is 
fixed  giving  the  increment  in  drag  as  depicted  by  the  preceding 
comparison.  Compare  also  the  drag  polar  for  the  13  %  airfoil 
under  transition  fixed  conditions  at  M  =  0.76  to  its  transition 
free  counterpart  in  Figure  33.  Some  contribution  to  the 
increased  drag  under  transition  fixed  conditions  must  come  from 
the  trip  drag  of  the  transition  strip. 

To  show  the  Reynolds  number  effect  on  drag  one  may 

compare  the  transition  fixed  data  at  Reynolds  number  of  Re/c  = 

6  — 

12.5  X  10  to  its  transition  fixed  counterpart  at  Re  =  6.7  X 

106  .  Drag  polars  for  the  13  %  and  21  %  thick  airfoils  in 

Figures  27  and  35  can  be  used  this  comparison.  Based  on  the  fact 

that  for  similar  boundary  layers  the  local  skin  friction  * 

1/Rp ( x )  7  ,  one  would  expect  lower  drag  values  for  the  high 

Reynolds  number  case.  This  is  certainly  borne  out  by  the 

comparison  of  the  transition  fixed  data  in  Figures  27  and  35. 

2.3.2  Drag  versus  Mach  number 

The  drag  polars  of  the  four  airfoils  are  also  cross 
plotted  in  the  form  of  wake  drag  CDW  as  a  function  of  the  Mach 
number  M  for  a  number  of  selected  cases.  Both  the  transition 
free  and  transition  fixed  data  at  Reynolds  numbers  of  6.7  X  10^ 
and  12.5  X  10^  are  plotted  at  Cr _  =  0.3,  0. 5,0.6  and  0.7.  The 
plots  are  presented  in  Figures  39  to  54. 
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Since  airfoil  drag  has  already  been  discussed  at  length 

in  the  previous  section  ,  only  some  interesting  features  of  drag 

£ 

at  Re/c  =  6.7  X  10  (  where  a  long  extent  of  laminar  flow  was 

observed)  will  be  highlighted  here.  Very  often  in  the  CLB  v  M 
type  of  plots  for  NLF  supercritical  airfoils  we  witness  the  drag 
'bucket  phenomenon',  which  as  the  name  suggests  represents  a 
noticeable  depression  in  the  drag  levels  close  to  the  design 
Mach  numbers.  Indeed  it  is  at  or  close  to  these  flow  conditions 
where  as  much  as  65  to  70%  laminar  flow  (10%  airfoil)  has  been 
observed  Ref. 7.  Also  see  section  2.5. 

For  comparison  of  drag  performances  between  airfoils  and 

g 

the  extent  of  the  bucket  region  for  the  case  Re/c  =  6.7  X  10  , 
consult  table  1  and  Figures  39  to  42. 


2.4  Drag  Comparison  Against  Other  Airfoils 


Drag,  at  design  conditions  at  the  lower  Reynolds 
numbers,  obtained  from  the  current  set  of  NLF  airfoils  as  seen 
in  Figure  55  is  considerably  less  than  all  the  2  D  airfoils 
tested  in  this  facility  to  date.  Hoerner's  turbulent  shock  free 
drag  curve  also  provided  in  Figure  55  represents  a  minimum 
envelope  for  all  of  the  other  airfoils  shown.  It  must  be 
realized  that  some  of  the  airfoils  included  in  this  comparison 


used  fixed  transition  strips  and  thus  lost  the  benefit  of  NLF, 
whereas  others  did  make  use  of  NLF  and  still  produced  drag 
levels  comparable  with  those  of  fixed  transition.  This  could 
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have  been  due  to  unfavorable  pressure  gradients  (  as  in  the 
'peaky'  type  airfoils  )  or  tunnel  turbulence  levels  at  the  time 
of  testing.  Recent  improvements  Ref.  11  to  the  NAE  5  ft  X  5  ft 
wind  tunnel  might  have  cleaned  up  some  of  this  turbulence.  The 
turbulence  levels  averaged  to  about  0.16%  at  the  conditions 
where  the  current  airfoil  tests  were  conducted. 


Inspection  of  Figure  55  shows  that  minimum  drag  levels 

given  by  our  four  airfoils  are  indeed  well  below  Hoerner's 

£ 

curve,  especially  at  Re/c  =  6.7  X  10  where  the  drag  levels  for 
the  four  airfoils  range  between  0.0057  and  0.0064. 

To  compare  with  still  other  data  we  consider  the  NACA 
66  series,  which  is  a  family  of  low  drag  airfoils.  It  must 
however  be  noted  that  the  NACA  66  series  airfoils  are  earlier 
designs  and  are  most  suited  for  low  speed  and  low  lift 
applications.  A  comparison  between  the  aerodynamic  performance 
of  the  present  NAE  airfoils  and  NACA  66  series  airfoils  is  shown 
in  table  2.  Instead  of  further  cluttering  up  Figure  55,  a 
minimum  drag  CDW  comparison  along  with  other  important 

aerodynamic  parameters  such  as  CL  (max)  and  3CL/3a  (at  design  M 
),  between  NAE  airfoils  and  the  airfoils  selected  from  NACA  66 
series  is  carried  out  in  table  2. 

2.5  FLOW  VISUALIZATION 


Some  flow  visualization  data  were  also  recorded  to 
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substantiate  the  extent  of  the  natural  laminar  flow.  The  flow 
visualization  was  obtained  by  spraying  the  airfoils  with  a  thin 
film  of  oil  containing  a  dye  which  fluoresced  in  ultraviolet 
light.  The  application  required  some  care  as  a  film  that  was  too 
thick  could  cause  premature  transition  which  would  alter  the 
airfoil  characteristics  markedly  at  low  Reynolds  numbers. 

Figure  56  shows  the  flow  visualization  record  of  the  16% 

f. 

thick  airfoil  at  Re/c  =  6.7  XlO  ,  M  =  0.75  and  CL  =  0.6  which 
correspond  to  the  conditions  where  drag  bucket  was  experienced 
from  the  wake  drag  data  (see  Fig  41)  .  Various  regions  of  the 
flow  are  described  on  the  figure.  It  appears  that  at  Re/c  =  6.7 
XlO  laminar  flow  extends  up  to  70  %  chord  which  coincides  with 
the  shock  location.  The  disturbances  due  to  the  pressure  taps 
caused  a  transition  wedge  of  sufficient  spanwise  extent  to 
influence  the  centerline  probe  of  the  wake  traverse  apparatus. 
Although  the  models  were  cleaned  very  thoroughly  between  runs, 
the  flow  visualization  record  showed  a  few  isolated  transition 
wedges  due  to  minute  particles  from  tunnel  contamination  during 
the  run. 

Figure  57  shows  a  flow  visualization  record  for  the  13% 
thick  airfoil  close  to  the  conditions  where  minimum  bucket  drag 
was  measured  for  this  airfoil  ,  see  Fig  41.  The  test  conditions 
corresponded  to  M  =  0.78,  Re/c  =  6.7  x  10^  and  <x  =0.33.  The 
two  black  strips  on  the  airfoil  were  the  tapes  used  to  cover  the 
pressure  taps  on  the  airfoil.  The  laminar  flow  in  the  wider 
clean  spanwise  portion  of  the  airfoil  runs  well  up  to  =65%  of 
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the  chord  length.  The  transition  wedge  emanating  from  the  corner 
of  the  tape  is  fairly  noticeable  up  to  this  station  (x/c  =  0.65) 
where  it  gets  muddled  into  the  overall  turbulent  flow.  Note  the 
presence  of  two  distinct  specks  of  dirt  in  the  segment  between 
the  two  strips  giving  rise  to  earlier  transition  wedges. 

For  the  10%  thick  airfoil,  the  flow  visualization  record 

6 

is  shown  in  Fig  58,  at  conditions  M  =  0.8,  Re/c  =  6.7  X  10  and 
a  =  1.39,  which  again  correspond  to  the  bucket  conditions  in 
Fig  41.  A  laminar  flow  length  of  well  up  to  60%  is  measured  in 
this  case.  Note  the  presence  of  the  transition  wedge  on  the 
right  hand  side  caused  by  a  dirt  particle. 

2.6  PRESSURE  DISTRIBUTIONS 

Figures  59a  and  59b  respectively  compare  the  transition 

fixed  and  free  pressure  distributions  on  the  four  airfoils  at 

design  conditions,  i.e.  CL  =  0.6,  the  Reynolds  number  of  Re/c  = 

6 

12.5  X  10  and  their  respective  design  Mach  numbers.  To  obtain 
long  lengths  of  laminar  flow  one  would  expect  the  favorable  or 
flat  roof  top  pressures  to  drive  the  shock  as  far  downstream  as 
possible.  This  is  best  represented  by  the  pressure  distributions 
on  the  upper  surface  of  the  10  and  13  %  airfoils.  The  10  % 
airfoil  was  able  to  delay  shock  until  66%  of  the  chord  under 
both  free  and  fixed  transition  conditions.  For  the  13  %  airfoil 
the  shock  location  was  about  50  %  under  the  design  transition 
fixed  conditions  and  about  45  %  under  off  design  transition  free 
conditions.  The  21  and  the  16  %  produced  favorable  pressure 


20 


distributions  only  up  to  about  2C  %  of  the  chord  followed  by  a 
gentle  adverse  pressure  gradient  (the  gradient  is  more 
noticeable  for  the  16  %  airfoil)  under  transition  fixed 
conditions.  Under  transition  free  conditions  the  21  %  airfoil 
seems  to  show  a  flatter  and  longer  roof  top  pressure 
distribution  when  compared  to  its  transition  fixed  counterpart. 
Note  also  a  more  pronounced  pressure  recovery  after  shock  for 
the  21%  airfoil  under  transition  free  conditions. 

The  pressure  levels  also  tend  to  be  a  little  higher 
under  fixed  transition  conditions  giving  rise  to  a  stronger 
shock  than  the  smearing  of  the  pressure  levels  for  the  free 
transition  case  producing  weaker  shock. 

The  bottom  surface  shows  favorable  pressure  gradient  for  all 
airfoils  up  to  a  distance  of  about  45  to  50  %  of  the  chord  for 
both  cases  considered  above. 

To  provide  another  comparison,  Figure  60  shows  the  free 
transition  data  at  Re/c  =  6.7  X  10®,  with  the  CL  and  Mach  number 
again  corresponding  to  the  design  conditions  of  the  airfoils. 

The  10  %  airfoil  still  continues  to  show  a  remarkably  flat  roof 
top  pressure  distribution  on  the  top  surface,  stretching  up  to 
about  65  %  of  the  chord  length  even  at  this  off  design 
transition  f-.ee  condition.  The  13  %  and  21%  airfoils  too 
maintain  a  steady  flat  pressure  distribution  up  to  about  45  %  of 
the  chord  length  indicating  the  presence  of  substantial  lengths 
of  laminar  flow.  The  16  %  airfoil,  however  shows  a  noticeable 
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adverse  pressure  gradient  past  the  20%  chord  length  signaling 
the  onset  of  relatively  earlier  transition.  The  bottom  surfaces 
have  once  again  maintained  favorable  pressure  gradients  up  to  45 
to  50  %  of  the  chord  lengths. 

Notice  the  presence  of  a  small  'kink'  in  the  16  %  pressure 
data  at  about  3  %  of  the  chord  length  on  all  three  figures  which 
was  very  likely  caused  by  the  presence  of  a  leaky  pressure  port. 

3.  COMPARISONS  OF  GRUMFOIL  WITH  THE  EXPERIMENTAL  DATA 


In  reference  6  some  comparisons  were  made  between  our 
experimental  data  on  the  present  airfoils  and  the  computer  codes 
BGK  f 12  ]  ,  GRUMFOIL  [13-14]  and  DRELA  [15].  The  comparisons 
indicated  that  GRUMFOIL  and  DRELA  were  quite  close  to  good 
agreement  with  GRUMFOIL  being  marginally  better.  Thus  in  the 
paper  we  will  consider  only  this  code  and  concentrate  our 
attention  on  the  10%  and  13%  airfoils.  We  will  also  include 
studies  on  the  free  transition  cases  as  well  as  the  fixed 
transition  cases  given  in  Reference  6.  It  was  found  that  in  the 
free  transition  cases  it  was  more  difficult  to  obtain  a  solution 
-  the  code  is  more  robust  if  transition  is  fixed  and  even  here  a 
good  initial  guess  for  the  solution  is  usually  needed.  To  this 
end  NAE  altered  the  GRUMFOIL  input  logic  to  accept  successive 
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starts  with  the  previous  final  solution.  For  the  first  input 
line  it  is  best  to  start  with  a  low  Mach  number  and  lift,  and 


transition  far  forward. 
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Thu  GRUMFOI L  code  has  3  options  for  determining 

transition.  These  are  based  on  Crabtree's  correlation  [16], 

Michel's  [17]  and  that  of  Stevens,  Gordia  and  Braden  [18].  The 

latter  two,  for  our  test  data,  did  not  give  convergent  solutions 

so  we  were  limited  to  using  Crabtree's  correlation.  Using 

experimental  data  at  Reynolds  numbers  up  to  7  million  based  on  3 

different  airfoils  and  a  flat  plate,  Crabtree  was  able  to  show 

that  a  reasonable  curve  through  the  mean  data  could  be  obtained 

by  plotting,  at  transition,  the  boundary  layer  Reynolds  number 

2 

U0/v  against  a  pressure  g. adient  parameter  U'9  /v.  The  GRUMFOIL 
code  uses  a  similar  curve  to  predict  transition. 

The  comparison  of  theory  with  experiment  must  be  done 
carefully  as  the  wind  tunnel  >‘>oh  number  cannot  be  defined  with 
100%  accuracy.  The  top  and  bottom  wall  correction  due  to  Mokry 
and  Ohman  [19]  has  been  used  at  NAE  for  several  years  and  is 
similar  to  corrections  applied  in  other  porous  wall  tunnels.  At 
NAE  we  have  always  assumed  that  sidewall  interference  is 
negligible  since  sidewall  suction  is  applied  in  the  vicinity  of 
the  airfoil.  However  recent  work  by  Lynch  [20-21]  indicates  that 
there  may  be  a  correction  due  to  sidewall  boundary  layer  as 
applying  the  suction  does  not  result  in  a  perfectly  parallel 
flow  at  the  walls.  Murthy  has  computed  a  correction  due  to 
sidewall  boundary  layer  variation  96  /9x  for  the  NAE  facility. 
His  work  [221  is  based  on  an  assumption  of  2-D  growth  of  the 
boundary  layer  with  no  lateral  variation.  This  assumption  is 
questionable  but  it  does  lead  to  a  correction  of  a  right  order 
of  magnitude  when  applied  to  test  data  by  Lynch  [20-21]  and 
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Jones  et  al  [6].  However  even  these  comparisons  are  questionable 
since  a  comparison  cannot  be  made  as  the  scaling  used  for 
is  in  error  (see  [6]).  In  the  comparisons  in  the  present  report 
we  will  always  use  the  ratio  p/pQ  since  these  are  wind  tunnel 
measured  quantities  which  are  invariant  with  free  stream 
conditions  M  and  p_ . 

The  difficulties  in  getting  a  good  match  are  shown  in 
Figure.  61.  Here,  transition  is  fixed  at  7%  and  15%  on  the  upper 
and  lower  surfaces  respectively.  It  can  be  seen  that  applying 
the  full  correction  indicated  by  Murthy  (AM  =  -0.015),  in  Fig. 
61a,  the  upper  surface  distribution  is  well  predicted  but  the 
lower  surface  pressures  are  too  low.  By  increasing  the 
computational  Mach  number  such  that  AM  =  -0.008  we  see,  from 
Fig.  61b,  that  the  lower  surface  matches  well  but  the  upper 
surface  is  not  well  predicted. 

An  interesting  feature  shown  in  Fig.  61  are  the 
'ripples'  in  the  upper  surface  pressure  distribution.  These  are 
thought  to  be  due  to  oblique  waves  emanating  from  the 
interference  of  the  sidewall  with  the  airfoil  leading  edge. 

These  have  been  noticed  by  Boeing  [23]  who  also  have  flow  viz 
pictures  showing  this  effect.  To  verify  that  it  was  not  a  'poor 
model'  effect  in  our  case  we  ran  a  omputation  using  design 
coordinates  and  measured  coordinates.  The  results  were  almost 
identical  indicating  that  the  model  was  accurately  made. 


Until  the  sidewall  effect  is  understood  more  clearly 
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( NAE  may  use  a  full  5'  span  model  to  minimize  sidewall  effects) 
we  will  use  in  this  report  the  correction  as  evaluated  according 
to  formulae  by  Murthy  [22]  and  shown  in  Fig.  62.  Some  typical 
surface  pressure  comparisons  are  shown  in  Fig.  63-65.  The  first 
set,  Fig.  63,  is  for  the  10%  foil  with  transition  fixed  at  7% 
and  15%  while  the  remaining  set  are  for  the  10%  (free 
transition)  and  13%  (fixed).  For  one  of  the  cases,  Fig  64c,  the 
convergence  and  boundary  layer  parameters  are  shown  in  Fig  66. 

Note  that  for  the  free  transition  case  on  the  13% 
airfoil  solutions  are  not  available.  Several  attempts  were  made 
to  run  GRUMFOIL  in  a  mode  where  transition  was  either  determined 
using  Crabtree's  criterion  or  fixed  at  the  specified  input 
value-  whichever  came  first  ( KT  =1,  XTANU  >  0,  X-TRANL  >0).  Since 
convergence  is  difficult  when  transition  is  far  back,  the  code 
was  operated  by  moving  the  transition  back  in  successive 
solutions,  while  this  approach  worked  in  some  cases  for  the  10% 
airfoil,  it  failed  in  our  attempts  for  the  13%.  It  seems  that 
during  the  iterations  the  transition  point  moved  backwards  and 
forwards  without  stabilizing.  It  may  be  that  a  more  gradual 
change  from  laminar  to  turbulent  flow  would  help  rather  than 
having  a  sudden  change. 

Generally  the  trends  shown  by  the  predictions  are  quite 
good  especially  in  the  fixed  transition  cases.  The  free 
transition  results  are  more  difficult  to  obtain  and  the  code 
will  often  not  converge  to  a  steady  state.  Thus  it  is  not 
possible  to  do  all  the  comparisons  one  would  like  in  this  case. 
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Even  so  in  the  cases  that  did  converge  most  of  the  trends  seem 
reasonable  (e.g.  transition  moving  downstream  as  the  pressure 
gradient  becomes  more  favorable)  and  the  code  shows  a  prediction 
of  the  'drag  bucket'  even  though  the  drag  level  is  too  high 
(Fig.  67).  This  drag  data  from  GRUMFOIL  was  obtained  using  R  = 

v 

6.7  million  with  about  0.6  and  increasing  the  Mach  number  in 

steps  from  0.489  to  0.802  (Mc  from  0.5  to  0.82).  The  computed 
transition  points  and  other  data  are  given  in  Table  3.  Also 
shown  in  Fig.  67  are  the  data  for  R  =11  million  and  16.7 
million.  It  can  be  seen  here  that  the  drag  trend  with  Reynolds 
number  is  not  correct  since  the  lower  Reynolds  number  data,  with 
transition  further  aft,  should  be  predicting  lower  drag.  Note 
that  the  data  is  not  complete  for  the  higher  Reynolds  numbers  as 
convergence  problems  are  encountered  at  these  conditions. 

4.0  CONCLUSIONS 

Four  NLF  airfoils  having  thicknesses  of  0.1,  0.13,  0.16  and  0.21 
have  been  investigated  both  theoretically  and  experimentally  in 
the  High  Speed  Aerodynamics  Laboratory  of  NAE. 

At  or  close  to  the  design  flow  conditions  as  much  as  70%  natural 
laminar  flow  has  been  observed  on  some  airfoils. 

The  minimum  drag  levels  (C^)  for  the  four  airfoils  range  from  a 
low  of  about  0.0045  to  a  high  of  about  0.0051. 

A  complete  comparison  of  experimental  results  with  GRUMFOIL 
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calculations  was  not  possible  due  to  non  convergence  of  the  code 
in  a  large  number  of  transition  free  cases.  The  cases  with  fixed 
transition  were  less  likely  to  fail.  The  available  results  for 
the  10%  airfoil  in  free  transition  indicate  a  correct  prediction 
of  the  drag  bucket  at  a  Mach  number  of  about  0.8. 

Having  investigated  the  NLF  concept  in  2  D,  there  is  a  need  to 
verify  some  of  these  results  in  3  D  by  similarly  investigating  a 
generic  reflection  plane  (preferably)  or  complete  model  in  NAE's 
5  ft  High  Reynolds  number  Trisonic  Facility.  Ultimately,  we  hope 
to  apply  this  NLF  c^n.^pt  to  a  real  flight  experiment  which  may 
involve  gloving  a  portion  of  the  wing  of  a  Flight  Research  T33 
trainer  aircraft. 
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TABLE  1:  SOME  INTERESTING  DRAG  FEATURES  OF  THE  FOUR  NLF 
AIRFOILS  AT  Re  =  6.7  x  106  ft"1  (TRANSITION  FREE) 


|  AIRFOILS  | 

1 

1 

1 

1 

|  CONDITIONS  |  CLB  =0.3 

lCLB  '  »• 

5ICLB  =  0 

•6I  CLB  = 

0-7  | 

NAE  68-060-21:1 

Bucket  | 0 . 66  <  M< 

| 0 . 66 <M  < 

|  0 . 6  6  <  M 

< | 0 . 64  <M 

< 

region  | 0 . 7 1 2 

| 0 . 712 

|  0 . 71 

|  0 . 71 

1 

CDW  ( min ) =  °-005;L  CL  a  0-497  at  M  =  0.68  and  Re  = 

6.7  X  106  ft-1 

Mdr  based  3CD/3M=0 . 1 | 0 . 703 

| 0 . 700 

|  0 . 712 

|0.689 

1 

CDW  (bucket  min)  | 0 . 0 0 5 7 

| 0 . 0058 

|0.0064 

|0.0072 

1 

NAE  72-060-16:1 

Bucket  | 0 . 7  <  M  < 

| 0 . 70 <M  < 

|  0 . 70<M 

<  | 0 . 72  <M 

<  1 

region  | 0 . 7  7 

|0.76 

|  0.76 

|  0 . 76 

1 

CDW  («i»)  '  °-°051  Cl  *  °-«4 

at  M  =  0. 

75  and  Re 

=  6.7  X  106 

ft-1 

Mdr  based  3CD/3M=0 . 1 | 0 . 764 

|0.752 

|0.748 

|0.758 

1 

CDW  (bucket  min)  |  0 . 0 0 5 7 

|0.0051 

|0.0056 

(0.0070 

1 

NAE  76-060-13:1 

Bucket  |  0 . 72  <  M  < 

| 0 . 74  <M  < 

|  0 . 7  2  <M 

<  | 0 . 72  <M 

<  1 

region  | 0 . 8  0 

|  0 . 80 

|  0 . 79 

|  0.80 

1 

CnM  ,  =  0.0045  Cr  =  0.300 

DW  (nun)  L 

at  M  =  0. 

76  and  Re 

=  6.7  X  106 

ft-1 

Mdr  based  3CD/3M=0 . 1 | 0 . 772 

|0.779 

(0.78 

|0.783 

1 

CQW  (bucket  min)  (0.0045 

10.0051 

(0.0056 

|0.0078 

1 

NAE  80-060-10:1 

Bucket  |  0 . 7 4  <  M  < 

|  0 . 7  4  <M  < 

| 0 . 74  <M 

<  | 0 . 74  <M 

<  1 

region  | 0 . 8  2 

|  0 . 82 

|0.82 

|  0 . 81 

1 

CnM  .  =  0.0045  CT  =  0.540 

DW  (mm)  L 

at  M  =  0. 

81  and  Re 

=  6.7  X  106 

ft-1 

Mdr  based  3CD/3M=0 . 1 | 0 . 812 

|0.811 

|0.819 

|0.800 

1 

CQW  (bucket  min)  | 0 . 0 0 5 2 

(0.0048 

10.0058 

|0.0080 

1 

Individual  drag  rise  Mach  number  MDR  for  each  case  is 
appropriately  identified  on  the  CDW  V  M  plots. 
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TABLE  2:  COMPARISON  WITH  THE  NACA  66  SERIES 


AIRFOILS 

t/c 

CL( des ) 

r 

L  (max) 

3CL/3a 

at  M  (des) 

r 

D  ( min 

NAE  68-060-21:1 

0.21 

0.60 

1.6800 

0.176 

0.0051 

NACA  66^  -221 

0.21 

— 

1.4900 

0.111 

0.0042 

NAE  72-060-16:1 

0.16 

0.60 

1.6400 

0.198 

0.0051 

NACA  662  -215 

0.15 

0.25 

1.4500 

0.106 

0.0035 

NAE  76-060-13:1 

0.13 

0.60 

1.6200 

0.196 

0.0045 

NACA  661  -212 

0.12 

0.15 

1.4500 

0.106 

0.0035 

NAE  80-060-10:1 

0.10 

0.60 

1.1900 

0.240 

0.0045 

NACA  66  -210 

0.10 

0.15 

1.2800 

0.110 

0.0033 
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TABLE  3:  GRUMFOIL  AND  EXPERIMENTAL  RESULTS  FOR  THE 
10%  AIRFOIL:  TRANSITION  POINTS  AND  DRAG  VALUES. 

Rc  =  6.7  x  106,  CL=  0.6 


M 

c 

C°w 

S 

CDF 

CD 

OO 

CD 

wave 

*r 

u 

\  (%) 

L 

0.5 

0.0075 

0.0076 

0.0051 

0.0071 

0.0005 

55 

0.7 

0.0074 

0.0087 

0.0055 

0.0081 

0.0006 

6 

54 

0.74 

0.0077 

0.0084 

0.0051 

0.0078 

0.0006 

12 

54 

0.76 

0.0067 

0.0081 

0.0049 

0.0072 

0.0009 

21 

54 

0.78 

0.0070 

0.0078 

0.0046 

0.0066 

0.0012 

36 

54 

0.79 

0.0060 

0.0078 

0 . 0044 

0.0064 

0.0014 

42 

54 

0.8 

0.0060 

0.0077 

0.0042 

0.0061 

0.0016 

48 

52 

0.81 

0.0057 

0.0067 

0.0034 

0.0058 

0.0009 

70 

52 

0.82 

0.0100 

0.0114 

0.0036 

0.0055 

0.0059 

70 

50 

NOTATION 

Cn  EXPERIMENTAL  WAKE  DRAG 

Dw 

GRUMFOIL: 


Cn  TOTAL  DRAG  INCLUD.  !G  BASE  DRAG 

c  SKIN  FRICTION 

°DF 

r  BASED  ON  MOMENTUM  THICKNESS  FOR  DOWNSTREAM 

Da 

C  DRAG  DUE  TO  THE  SHOCK  WAVE 

uwave 

x  TRANSITION  POINTS  ON  THE  UPPER  AND  LOWER 

T  SURFACES 


to 


POTENTIAL  LAMINAR  RUN  ON  UPPER  SURFACE (FREE 
**  H  TRANSITION) 

TDAucTTtnM  »  CHORDWISE 

TRANSITION  ^  ,  EXTENT 

/  _  -  "  "  »  DEPENDS  ON 

- T~A  l/t/c 

mdes  W 

/^FORtti^ 


POTENTIAL  LAMINAR  RUN 

M - — - 

LOWER  SURFACE  I 


FIG.  2:  FEATURES  OF  DESIGN  VELOCITIES 
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FIG.  4:  COEFFICIENT  OF  LIFT  VERSUS  ANGLE  OF  ATTACK 

C,  n  VERSUS  a 


060-21 :1 
060-16:1 
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VERSUS  M,  Re  =  6.7  x  10b  (FREE  TRANSITION) 


A  NAE  72-060-16:1 


INEAR  COMPRESSIBILITY  VERSUS  , - -  (FREE  TRANSITION) 


VERSUS  M,  Re  =  6.7  x  106  (TRANSITION  FIXED) 


060-21:1 

060-16:1 

060-13:1 

060-10:1 
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VERSUS  M,  Re  =  12.5  x  106  (TRANSITION  FIXED) 


S  Zl  t-U  90l. 


LINEAR  COMPRESSIBILITY  *3—  VERSUS  ~j=  (TRANSITION  FIXED) 


-  THICKNESS  OF  AIRFOIL 
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FIG.  12:  MAXIMUM  LIFT  AND  LIFT  AT  SEPARATION  VERSUS  MACH  NUMBER 
C,  ,  C,  VERSUS  M,  Re  =  12.5  x  106  (FREE  TRANSITION) 

‘-MAX’  USEP 


-  THICKNESS  OF  AIRFOIL 
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MAXIMUM  LIFT  AND  LIFT  AT  SEPARATION  VERSUS  MACH  NUMBER 
L  ,  C.  PD  VERSUS  M,  Re  =  6.7  x  106  (TRANSITION  FIXED) 


-  THICKNESS  OF  AIRFOIL 


48 


FIG.  14:  MAXIMUM  LIFT  AND  LIFT  AT  SEPARATION  VERSUS  MACH  NUMBER 
Cl  .  CU  VERSUS  M,  Re  =  12.5  x  106  (TRANSITION  FIXED) 

MAX  SEP 


-0.08 
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Re  =  6.7  x  106,  CL_  =  0.6,  (FREE  TRANSITION) 
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FIG.  16:  PITCHING  MOMENT  CM  VERSUS  MACH  NUMBER  M, 
Re  =  6.7  x  106,  CLr  =  0.6,  (TRANSITION  FIXED) 
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FIG.  19:  THE  WAKE  RAKE 
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FIG.  20:  WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CDw  VERSUS  CLb,  (FREE  TRANSITION) 

Re  =  6.7  x  106,  M  =  0.2  &  0.3 


WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CLw  VERSUS  CLq,  (FREE  TRANSITION) 

Re  =  6.7  x  106,  M  =  0.6 


FIG.  23:  WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CDw  VERSUS  CLb  (FREE  TRANSITION) 

Re  =  6.7  x  106,  m  =  0.7 


0.032 


58 


FIG.  24:  WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CDw  VERSUS  CLb,  (FREE  TRANSITION) 

Re  =  6.7  x  106,  M  =  0.76,  0.74 


FIG.  25:  WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CDw  VERSUS  CLb,  (FREE  TRANSITION) 

Re  =  6.7  x  106,  M  =  0.80,  0.81 
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FIG.  26:  WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CDw  VERSUS  CLb  (TRANSITION  FIXED) 

Re  =  6.7  x  106,  M  =  0.2 


0.032 
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FIG.  28:  WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CDw  VERSUS  CLb  (TRANSITION  FIXED) 

Re  =  6.7  x  106,  M  =  o.68  -  0.72 
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FIG.  30:  WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CDw  VERSUS  CLb,  (FREE  TRANSITION) 

Re  =  12.5  x  106,  M  =  0.5 


0.032 


WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CDw  VERSUS  CLb,  (FREE  TRANSITION) 

Re  =  12.5  x  106,  M  =  0.6 


FIG.  34:  WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CDw  VERSUS  CLb,  (FREE  TRANSITION) 

Re  =  12.5  x  106,  M  =  0.78,  0.80 


WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CDw  VERSUS  CLb>  (TRANSITION  FIXED) 

Re  =  12.5  x  106,  m  =  0.6 


FIG.  36:  WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CDw  VERSUS  CLb,  (TRANSITION  FIXED) 

Rs  =  12.5  x  108,  M  =  0.70  -  0.72 


WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CDw  VERSUS  CL_,  (TRANSITION  FIXED) 


FIG.  38:  WAKE  DRAG  COEFFICIENT  VERSUS  LIFT  COEFFICIENT 
CDw  VERSUS  CLb,  (TRANSITION  FIXED) 


0.0125 


FIG.  39:  WAKE  DRAG  COEFFICIENT  VERSUS  MACH  NUMBER 
CDw  VERSUS  M,  (FREE  TRANSITION) 

CLb  =  0.3,  Re  =  6.7  x  106 
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VERSUS  M,  (FREE  TRANSITION) 


0.0125 


FIG.  41:  WAKE  DRAG  COEFFICIENT  VERSUS  MACH  NUMBER 
CDw  VERSUS  M,  (FREE  TRANSITION) 

CLd  =  0.6,  Re  =  6.7  x  1Q6 


0.0125 


FIG.  42:  WAKE  DRAG  COEFFICIENT  VERSUS  MACH  NUMBER 
CDw  VERSUS  M,  (FREE  TRANSITION) 
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WAKE  DRAG  COEFFICIENT  VERSUS  MACH  NUMBER 
CDw  VERSUS  M,  (TRANSITION  FIXED) 


NAE  68-060-21:1 
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VERSUS  M,  (FREE  TRANSITION) 


0.024 


F!G 


MACH  NUMBER  SHOWN  AS  INTEGER 

ALL  AIRFOILS  TESTED  AT  NAE  5  FT  X  5  FT  FACILITY 


FIG.  55:  DRAG  COMPARISONS  WITH  OTHER  AIRFOILS  USING  CD 
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1  TRANSITION  WEDGE  DUE  TO  LINE  OF 

PRESSURE  TAPS 

2  TRANSITION  WEDGES  DUE  TO  SPECKS 

OF  DIRT 

3  LAMINAR  RUN  BACK  TO  70%  CHORD 

4  TURBULENT  FLOW  IN  PRESSURE 

RECOVERY  REGION 

5  SHOCK  LOCATION  IN  TURBULENT  FLOW 

6  PRESSURE  TAPPINGS  OFFSET  TO 

REDUCE  EFFECTS  OF  TRANSITION 
ON  PROBE  NO.  1 


FLOW  VISUALIZATION  OF  NAE  72-060-16:1  AT  Re  =  8  MILLION 
M  =  0.75,  AND  CL  =  0.6 


FIG.  56: 


FIG.  57:  FLOW  VISUALIZATION  OF  NAE  78-060-13:1 
AT  M  =  0.78,  Re  =  6.7  x  106,  a  =  0.33 
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FIG.  58:  FLOW  VISUALIZATION  OF  NAE  80-060-10:1 
AT  M  =  0.8,  Re  =  6.7  x  106,  a  =  1.39 


LAMINAR 


FIG.  59a:  PRESSURE  DISTRIBUTION  AT  THE  DESIGN  CONDITIONS, 


9' 


FIG.  59b:  PRESSURE  DISTRIBUTION  UNDER  FREE  TRANSITION  CONDITIONS 

Re/c  =  12.5  x  106 
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FIG.  60:  PRESSURE  DISTRIBUTION  UNDER  FREE  TRANSITION  CONDITIONS, 

Re/c  =  6.7  x  106 
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RUN  31313  SECN  2  SCAN  4 
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0.0  0.2  0.4  0.6  0.8  1.0 
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FIG.  61a:  GRUMFOIL  VERSUS  EXPERIMENT.  'SIDEWALL' 
CORRECTION  OF  -0.015  APPLIED.  TRANSITION  FIXED  10%  AIRFOIL 
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FIG.  61b:  GRUMFOIL  VERSUS  EXPERIMENT.  'SIDEWALL* 
CORRECTION  OF  -0.008  APPLIED.  TRANSITION  FIXED  10%  AIRFOIL 
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FIG.  63a:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSTION  FIXED.  10%  AIRFOIL.  CL  »  0.298 
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RUN  31313  SECN  2  SCRN  3 
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FIG.  63b:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSITION  FIXED.  10%  AIRFOIL.  CL  =  0.471 
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RUN  31313  SECN  2  SCAN  5 
M  CLP  RN  RLP  CMP  COT  MMRX  MSH  XSH 
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X 

FIG.  63c:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSITION  FIXED.  10%  AIRFOIL.  CL  -  0.663 
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FIG.  63d:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSITION  FIXED.  10%  AIRFOIL  CL  =  0.796 
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FIG.  64a:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSITION  FREE.  10%  AIRFOIL.  CL  -  0.333 
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FIG.  64b:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSITION  FREE.  10%  AIRFOIL.  CL  »  0.486 
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RUN  31234  SECN  1  SCfiN  4 
M  CLP  RN  ALP  CMP  COT  MMflX  MSH  XSH 
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FIG.  64c:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSITION  FREE.  10%  AIRFOIL.  CL  =  0.601 
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M 

CLP  RN 

RLP  CMP 

COT 

MMRX 

MSH 

XSH 

0.799 

0.697  6.8 

1 .23  -0.146 

0.0076 

1.24 

1  .24 

0.69  *  * 

CXP 

0.783 

0.715  6.8 

-0.11  -0.172 

0.0075 

1.24 

1.24 

0.77  - 

GRM 

0.0  0.2  0.4  0.6  0.8  1.0 

X 

FIG.  64d:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRTANSITION  FREE.  10%  AIRFOIL.  CL  »  0.697 
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FIG.  64e:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSITION  FREE.  10%  AIRFOIL.  CL  «  0.767 
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FIG.  64f:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSITION  FREE.  10%  AIRFOIL.  CL  »  0.810 


P/PO 

0.9  0.8  0.7  0.6  0.5  0.4  0.3  0.2 


109 


RUN  31281  SECN  1  SCRN  5 


M 

CLP 

RN  RLP  CMP 

CDT 

MMRX 

MSH 

XSH 

0.762 

0.376 

12.7  -0.16  -0.139 

0.0100 

1.10 

0.00 

0.00  *  * 

EXP 

0.747 

0.384 

12.5  -0.29  -0.134 

0.0106 

1.06 

0.00 

0.00  - 

GRM 

0.0  0.2  0.4  0.6  0.8  1.0 


X 

FIG.  65a:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSITION  FIXED.  13%  AIRFOIL.  CL  «  0.376 
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FIG.  65b:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSITION  FIXED.  13%  AIRFOIL.  CL  «  0.522 


P/PO 

0.9  0.8  0.7  0.6  0.5  0.4  0.3  0.2 


111 


RUN  31281  SECN  1  SCAN  7 


n 

CLP 

RN 

RLP  CMP 

COT 

MMRX 

MSH 

XSH 

0.762 

0.618 

12.7 

0.99  -0.143 

0.0106 

1.22 

1.19 

0.53  *  ’ 

EXP 

0.747 

0.630 

12.5 

0.79  -0.133 

0.0125 

1.20 

1.18 

0.53  - 

GRM 

0.0  0.2  0.4  0.6  0.8  1.0 


X 

FIG.  65c:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSITION  FIXED.  13%  AIRFOIL  CL  *  0.618 
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FIG.  65d:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSITION  FIXED.  13%  AIRFOIL  CL  -  0.717 
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FIG.  65e:  COMPARISON  OF  GRUMFOIL  AND  EXPERIMENT, 
TRANSITION  FIXED.  13%  AIRFOIL.  CL  »  0.798 
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FIG.  66a:  CONVERGENCE  OF  GRUMFOIL  FOR  CASE  IN  FIG.  64c 
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FIG.  66b:  8*  COMPUTED  BY  GRUMFOIL  CODE  FOR  CASE  IN  FIG.  64c 
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FIG.  66c:  SHAPE  FACTOR  PREDICTED  BY  GRUMFOIL  FOR  CASE 

IN  FIG.  64c 
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FIG.  66d:  MOMENTUM  THICKNESS  PREDICTED  BY  GRUMFOIL 

FOR  CASE  IN  FIG.  64c 
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FIG.  66e:  SKIN  FRICTION  AS  PREDICTED  BY  GRUMFOIL  FOR  CASE 

IN  FIG.  64c 
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FIG.  66f:  ENTRAINMENT  COEFFICIENT  PREDICTED  BY  GRUMFOIL 

FOR  CASE  IN  FIG.  64c 
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FIG.  67:  A  COMPARISON  OF  EXPERIMENTAL  AND  GRUMFOIL  D 

VALUES  CL  =  0.6.  10%  AIRFOIL 
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